Net-Baryon Physics: Basic Mechanisms 
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It is well known that, in nuclear collisions, a sizable fraction of the available energy is carried 
away by baryons. As the baryon number is conserved, the net-baryon B — B retains information 
on the energy- momentum carried by the incoming nuclei. A simple but consistent model for net- 
baryon production in high energy hadron-hadron, hadron- nucleus and nucleus-nucleus collisions is 
presented. The basic ingredients of the model are valence string formation based on standard PDFs 
with QCD evolution and string fragmentation via the Schwinger mechanism. The results of the 
model are presented and compared with both data and existing models. These results show that a 
good description of the main features of net-baryon data is possible on the framework of a simplistic 
model, with the advantage of making the fundamental production mechanisms manifest. 

PACS numbers: 12.38.Aw, 12.39.-x, 12.40.Nn, 13.85.Ni, 24.85.+p 



I. INTRODUCTION 

In hadron-hadron, hadron-nucleus and nucleus-nucleus 
interactions a sizable fraction of the available energy in 
a collision is carried away by baryons [l|, H, Q • ^^ ^^^ 
baryon number is conserved, the measured net-baryon, 
B — B^ keeps track of the energy- momentum carried by 
the incoming particles. An important question to be 
asked is: how does the fraction of energy carried by the 
net-baryon evolve as a function of the centre-of-mass col- 
lisional energy per nucleon, y^ ? This question is im- 
portant because, on one hand, a decrease of the fraction 
of energy going into the net-baryon implies more energy 
available to create the deconfined quark-gluon state of 
matter, and on the other hand, such a decrease may re- 
duce the possibility of producing fast particles in very 
high energy cosmic ray experiments. For more than 30 
years, since the ISR at CERN, particle production stud- 
ies have been limited to mid rapidity. Fortunately, with 
RHIC, large rapidity data became available, and, hope- 
fully, the same will happen for the LHC. In fact, if one 
does not measure the physics at high rapidity the most 
elementary physical constraint, namely energy conserva- 
tion, cannot be applied 0, Q . 

In most of the existing Monte Carlo models [y, 0, H, Q , 
the physics of net-baryon production is very much ob- 
scured by the complexity of extensive and detailed codes. 
Often, the basic production mechanisms do not appear 
in a transparent way. In this paper we present a simple 
but consistent model for net-baryon production in high 
energy hadron-hadron (h-h), hadron-nucleus (h-A) and 
nucleus-nucleus (A-A) collisions. As it happens in most 
of the existing codes, we shall work in the framework of 
the Dual Parton Model (DPM) [T^ with string formation 



(valence strings, in the present case) based on standard 
Parton Distribution Functions (PDFs) with QCD evolu- 
tion and string fragmentation via the Schwinger mecha- 
nism. The basic ingredients of the model are: 

• Formation of extended color fields or strings, mak- 
ing use of PDFs for valence quarks; 

• Evolution with momentum transfer Q^, as a conse- 
quence of the QCD evolution of the PDF's; 

• Fragmentation of each string with formation of a 
fast baryon (net-baryon) and other particles. 

This paper is organized as follows. In section [III the pre- 
dictions of current models and the available experimental 
data are briefly reviewed and compared. In section HIH 
our simple net-baryon model is described. The results 
of the model are presented and compared with data and 
with existing models in section HVl A summary and brief 
conclusions are presented in section \V\ 



II. CURRENT MODEL PREDICTIONS AND 
EXPERIMENTAL DATA 

The data presently available on net-proton or net- 
baryon production is scarce. The most recent results are 
from RHIC [ll|. These results are presented both in 
terms of net-baryon rapidity y and in terms of rapidity 
loss, defined as 



i^y) = yp- iv) 



(1) 



where yp is the beam rapidity and {y) is the mean net- 
baryon rapidity after the collision, given by 
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FIG. 1: The net-proton distribution at AGS [13] (Au+Au at 
Vi ^ 5 GeV), SPS [3 (Pb+Pb at ^s - 17 GeV), and RHIC 
(Au+Au at y/s — 200 GeV) is shown (upper plot). The data 
correspond to the 5% most central collisions and the errors 
are both statistical and systematic (the light gray band shows 
the 10% overall normalization uncertainty on the E802 points, 
but not the 15% for E917). The data have been symmetrised. 
For RHIC data, filled points are measured and open points are 
symmetrised, while the opposite is true for AGS and SPS data 
(for clarity). At AGS weak decay corrections are neg ligible 
and at SPS they have been applied. Taken from [Tl|. The 
distribution of A = ^ — ^/p is shown (lower plot) for SPS and 
AGS data. 



Here, N^art is the number of participants in the collision 
and N^_^ is the net-baryon number. It is worth not- 
ing that the net-baryon results depend directly on the 
number of participants in the collision and thus on the 
collision centrality. The relation between the number of 
participants and the impact parameter of the collision 
can be established, for instance, in the context of the 
Glauber model [l2l |. 

The net-baryon rapidity distributions at different 
centre-of-mass energies for different beam-target systems 
were compared in [ll[ . In figure [H the net-proton dis- 
tributions at AGS (Au+Au at ^/s c^ 5 GeV) [l3, SPS 
(Pb+Pb at V^ - 17 GeV) [13, and RHIC (Au+Au at 
y/s = 200 GeV) [ll] for central collisions (the 5% most 
central) are shown (upper plot). The distributions show 
a strong energy dependence, with the mid rapidity region 
corresponding to a peak at AGS, a dip at the SPS, and a 
broad minimum at RHIC. The distribution of A = y — yp 



FIG. 2: Rapidity loss as a function of beam rapidity (in the 
CM). The hatched area indicates the unphysical region, and 
the dashed line shows the phenomenological scaling (Sy) = 
0.58yp. The inset plot shows the BRAHMS net-baryon dis- 
tribution (data points) with fits (lines) needed to extrapolate 
up to the beam rapidity, explaining the large uncertainty as- 
sociated to the BRAHMS rapidity loss measurement. Taken 
from [HI. 



is shown (lower plot) for SPS and AGS data. It is appar- 
ent that neither Feynman scaling nor limiting fragmen- 
tation are satisfied in net-baryon production. 

In figure O the rapidity loss data is shown as a func- 
tion of the beam rapidity (in the centre-of-mass frame) 
for AGS, SPS and RHIC. The large uncertainty associ- 
ated with the BRAHMS data point is, as illustrated in 
the inset plot, related to a relatively unconstrained ex- 
trapolation to the high rapidity region, required for the 
rapidity loss calculation. 




NA49 Data 
■ EPOS 1 .61 - witin weak decays 
EPOS 1 .61 - no weak decays 
OGSJET-II - no weak decays 



FIG. 3: The net-proton rapidity distributions from EPOS 
1.61 (full lines) and QGSJET-II.03 (dashed line) for Pb-Pb 
collisions at a/s ~ 17 GeV are shown and compared with data 
NA49 SPS data (points). For EPOS, the results obtained 
including (thick line) and excluding (thin line) the strangeness 
contribution are shown. For QGSJET, the "no-decays" curve 
is shown. For details on the data see figure [1] 



Let us now turn to net-baryon production as imple- 
mented in the existjjig Monte Carlo models. QGSJET- 
II Q and EPOS [8] are amongst the presently most 
widely used hadronic models in high energy and cosmic 
ray physics. To our knowledge, there is no systematic 
study comparing the predictions of these two models on 
net-baryon production between themselves or with ex- 
perimental data. 

In figure [31 net-proton rapidity distributions obtained 
with QGSJET-II.03 and EPOS 1.61 for Pb-Pb collisions 
at v^ :^ 17 GeV are shown and compared with exper- 
imental data. According to [ij] weak decay corrections 
have been applied to the data. Following [15], the impact 
parameter range (0 to 3.1 fm) corresponding to the cen- 
trality cut of 5% applied to data was selected in QGS JET 
and EPOS. For EPOS, the results obtained leaving all 
particles free to decay and contribute to the net-baryon 
are shown together with those obtained switching off all 
weak decays. The "all-decays" EPOS curve reproduces 
the trend seen in data, with an excess that could be due to 
the presence of net-baryon from weak decays. However, 
the strangeness effect compensation in the "no-decays" 
curve seems to be too strong. For QGSJET-II, the "no- 
decays" curve is shown (and was found to be very simi- 
lar to what is obtained from the model by default). One 
should recall that QGSJET-II is not expected to perform 
very well at such low energies. 



III. THE MODEL 
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FIG. 5: The two main valence string fragmentation diagrams. 



momentum is considered within this model. By choice, 
xi > X2 > xs, with xs < 0. The energy and momentum 
of each string are obtained adding directly the energy 
and momentum carried by each of the valence quarks. 
For string A 



Let us now describe a simple model for net-baryon pro- 
duction. In the spirit of [2, |3|, [lO, [iQ] , the basic assump- 
tion is that net-baryon production in proton-proton col- 
lisions is strongly correlated with the formation and frag- 
mentation of two color singlet strings, each one with two 
valence quarks from one of the protons, and one valence 
quark from the other proton. This is schematically shown 
in figure [H 





FIG. 4: Schematic representation of a proton-proton collision, 
with the formation of two valence strings. 



Referring to the figure, let xi, X2 and xs be the frac- 
tions of momentum carried by the valence quarks forming 
string A. Quarks 1 and 2 are from the proton with pos- 
itive momentum in the proton-proton reference frame, 
and quark 3 is from the other proton. No transverse 



Estring = (^1 + ^^2 + (-^s)) 
Pstring = (^1 + ^^2 - (-^s)) 



2 
2 



M string = \/(xi + X2) (-X3) S, 



(3) 

(4) 
(5) 



where the momentum fractions xi, X2 and x^ are deter- 
mined using the quark PDFs. For each ^/s^ an effective 
momentum transfer Q^ adjusted from experimental data 
was chosen (see section |IVl for details). In this work the 
CTEQ6M parton distribution functions [17] were used. 

The simplest possible model for fragmentation is as- 
sumed. Each string decays into a baryon and a meson 
in the following way: the string is cut in two pieces and 
a qq pair is formed, from the vacuum, either between 
quarks 2 and 3 (belonging the different protons) or be- 
tween quarks 1 and 2 (belonging to the same proton, with 
positive momentum in this string A example). The quark 
pair with the largest momentum difference is chosen. The 
string piece that inherits two valence quarks originates 
the baryon, whereas the string piece that inherits one va- 
lence quark originates the meson. This mechanism cor- 
responds to the diagrams represented in figure [5l The 
weights of the two diagrams are, in this simple model, 
given only by kinematics. For string A the first diagram 
will be more probable (especially for low ^/s). However, 



the weight of the second diagram can be as large as 40% 
above LHC energies. 

The qq pair formed from the vacuum was taken to be 
either a iiiZ or a dd^ and the fuh quark combinatorics 
was then performed in order to determine the nature of 
the possible outcoming baryon. Both fundamental and 
excited states were considered, taking spin-dependent 
weights (2j+l). The decays of the unstable baryons were 
then performed and the outcoming nucleons included in 
the net-baryon calculations. The contribution from s 
quarks was not considered here. It was estimated from 
the Schwinger model to be about 25%. 

The string mass and momentum distributions as func- 
tion of Q^ are given in figure [6l The mass distribution 
applies to both valence strings, while the momentum dis- 
tribution for the case of string B is simply symmetric. 
The net-baryon rapidity distributions for different values 
of Q^ are presented in figure [7J for two different centre- 
of-mass energies. This figure illustrates not only the evo- 
lution with Q^ but also the clear effect of kinematics, 
due to the increase of y^, on the main features of the 
distribution. 

It is worth noting that the inclusion of diagram 2, in 
addition to diagram 1, with weights determined simply 
by kinematics, reproduces some of the effects predicted 
in models with string junctions [18, 1^ or popcorn [2^ 
mechanisms for the transport of baryon number from the 
beam rapidity into the central region y ^ 0. These effects 
can thus be achieved in a simple DPM model, based on 
valence strings and a Q^ parameterisation. 

For the purpose of estimating the net-baryon produc- 
tion in A- A collisions, we shall use the simple approxi- 
mation that dn/dy{B — 5) at a given ^/s is proportional 
to the number of participating nucleons, N parti 
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dy 
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dn 
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(6) 



This relation is presumably reasonable as net-baryon 
production is originated from valence strings associated 
to wounded (or participating) nucleons. At mid rapid- 
ity, saturation effects originate additional production of 
mesons and baryons. This affects {B-\-B), but not B — B 
production. An attempt to estimate nuclear effects cor- 
rection factors for the valence quark PDFs was made us- 
ing EKS98 HH and nDS [22^. Values below 10-15% were 
found. 



IV. RESULTS 

In order to fully define the model, we need to choose 
an effective Q^ value as input for the quark PDFs. This 
is done by adjusting the results of the model to the ex- 
perimental data, i.e., by choosing, for each ^/s^ the effec- 
tive Q^ for which the model better describes the data. 
In the fitting procedure, the global normalisation factor 
(the number of participants, relying on the approach of 
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FIG. 6: Evolution with Q^ of the string mass (upper plot) 
and momentum (lower plot) distributions for ^/s = 200 GeV. 
The mass distribution applies to both valence strings, while 
the momentum distribution for the case of string B is simply 
symmetric. 



yi (GeV) Collision Q^ 



N, 



part 



5 Au-Au 0.36l°;j; 260 ± 20 

17 Pb-Pb OMt^oll 300 ±20 

200 Au-Au 1.69tg;^8 280 ± 20 

TABLE I: Results of the fit to the effective Q^ and the number 
of participants at the different centre-of-mass energies. 



eq. dni is left as a second free parameter. In order to 
compare the results of this model on net-baryon to the 
experimental data on net-proton, we consider that net- 
proton is roughly 1/2 of {B — B). 

Figure [8] illustrates, for a specific y^ and fixing the 
number of participants, the ability of the model to gener- 
ate very different net-proton rapidity distributions, sim- 
ply by choosing different values of the effective Q^. At 
V^ = 200 GeV and ^/s c^ 17 GeV all the data points (see 
figure [1]) were included in the fit. At y^ c:^ 5 GeV, and 
due to energy conservation constraints imposed in the 
model, only the points up to the nominal beam rapidity 
were considered. 

The values obtained for the effective Q^ and the num- 
ber of participants Npart at the different centre-of-mass 
energies are given in table H 



values obtained for the number of participants. 





FIG. 7: Evolution of the net-baryon rapidity with Q^ at a/s -■ 
17 GeV (upper plot) and a/s = 200 GeV (lower plot). 
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FIG. 8: Net-baryon predictions of the present model for differ- 
ent values of Q^, and comparison with data, at a/s ^ 17 GeV. 



It is worth noting that the fitted N^art values are roughly 
compatible with what is expected from the literature, al- 
though about 20% lower. In fact, the average number of 
participants for Pb-Pb collisions at SPS energies is given 
in [131 as 362. Concerning Au-Au collisions, the value is 
expected to be only slightly lower (estimated as 344 to 
357 at RHIC energies [lil, S, Q) and expected to de- 
pend weakly on y^, in fair agreement with the obtained 
results. As pointed out in section lTTTl in this simple model 
strangeness contribution to the net-baryon was not con- 
sidered, and it amounts to about 25%, accounting for the 
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FIG. 9: The results of the present model for net-baryon rapid- 
ity (upper plot) and rapidity loss (lower plot) are compared 
to experimental data at different centre- of- mass energies. See 
figure [1] for details on the data points. 



Figure [9] shows the results of the present model for net- 
baryon rapidity distribution and rapidity loss, in compar- 
ison with experimental data at different centre-of-mass 
energies. A fair agreement is found at the measured en- 
ergies. For the lowest y^, the points above the nominal 
beam rapidity were not considered in the fit. 

We can now try to find a relation between the effec- 
tive Q^ and ^Js. The effective Q^ corresponds to the 
typical transverse size (area) of the parton (here, the va- 
lence quark). It is reasonable to assume, as in Regge 
phenomenology [25], that the average number of partons 
in a nucleon increases as a power of the centre of mass 
energy ^/s. Thus, ^fsjQ^ ~ i?^, where Rh is the nucleon 
radius which we take as fixed. It then follows that Q^ 
should grow according to 



Q' = Qo 



Vi 



A^ 



\GeV^ 



(7) 



The exponent A^ was determined by fitting eq. ([7j) to 
the available data points in table H 

The results are shown in figure [TOl where the data 
points are the Q^ values adjusted above, the line is the 
fit with eq. ([7j) and the shaded areas correspond to la 
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FIG. 10: The effective Q^ values at different centre-of-mass 
energies chosen by tuning the model to the experimental data 
are shown. The line shows the fit to the points using eq. (|7| 
and the shaded areas correspond to la variations of the fit 
parameters. 



FIG. 12: Evolution of the fraction of energy carried by the 
Net-Baryon with a/s. The prediction of the present model 
(with a shaded band corresponding to the la variation of the 
fit parameters) is shown, together with the results obtained 
with EPOS and QGSJET-II. The data point corresponds to 
the RHIC estimate given in |23| . 



variations of the fit parameters. The obtained results 
are A^ = 0.42to;?20 and Ql = OM^omt taking y^ = 
5 GeV. It should however be noted that we have just 
three points in the fit and the la band is rather wide. 
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FIG. 11: Predictions of the present model for net-baryon ra- 
pidity (upper plot) and rapidity loss (lower plot) including 
higher energies. 



Taking the obtained value of A^, the predictions of 
the present model for net-baryon rapidity and rapidity 
loss at higher centre-of-mass energies were obtained and 



are shown in figure [TTl covering both the LHC and the 
high energy cosmic ray regions. In addition, the fraction 
of the centre-of-mass energy carried by the net-baryon 
as a function of y^ was computed and is shown in fig- 
ure [12] (together with the Icr bounds). The predictions of 
EPOS 1.61 and QGSJET-II.03 are also shown. Accord- 
ing to [23], RHIC data point to about 27% of the initial 
energy remaining in the net-baryon after the collision. 
This result is also shown in the figure. It is apparent that 
EPOS, QGSJET-II and the present model show rather 
different trends. The present model and EPOS agree 
reasonably well at RHIC energies, where they are both 
slightly below the measured value. At these energies, 
QGSJET-II is somewhat above, but equally compatible 
with data. At higher energies, while in EPOS 1.61 the 
net-baryon is essentially zero, in the present model and 
in QGSJET-II.03 a sizable amount of energy is still as- 
sociated to the net-baryon. It should be noted that high 
energy effects such as string percolation may change these 
predictions |26l] . 



V. SUMMARY AND CONCLUSIONS 

A simple but consistent model for net-baryon produc- 
tion in high energy hadron-hadron, hadron-nucleus and 
nucleus-nucleus collisions was presented. The basic in- 
gredients of the model are valence string formation based 
on standard PDFs with QCD evolution and string frag- 
mentation via the Schwinger mechanism. The results of 
the model were presented and compared with data from 
RHIC, SPS and AGS, and with existing models, namely 
QGSJET-H and EPOS 1.61. The obtained results show 
that a good description of the main features of net-baryon 
data is achieved on the basis of this simple model, in 
which the fundamental production mechanisms appear 
in a transparent way. 



The free parameters in the model, the effective Q^ and 
the number of participating nucleons, were fitted to net- 
baryon data at y/s c^ 5, 17 and 200 GeV. A good fit 
to data is obtained within this model, and the values 
obtained for the number of participants are in agreement 
with what we would expect from the literature. 

Using the (scarce) net-baryon data at different centre- 
of-mass energies, a relation between the effective momen- 
tum transfer and the centre-of-mass energy was moti- 
vated and a prediction was obtained and extrapolated to 
higher energies for the evolution with ^/s of the fraction 
of the initial energy carried away by the net-baryon. A 
sizable amount of energy may be associated to the net- 
baryon, even at high energies. 
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